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Abstract 

Nonlinear evolution of circularly polarized Alfven waves are discussed by us- 
ing the recently developed Vlasov-MHD code, which is a generalized Landau- 
fluid model. The numerical results indicate that as far as the nonlinearity in 
the system is not so large, the Vlasov-MHD model can validly solve time evo- 
lution of the Alfvenic turbulence both in the linear and nonlinear stages. The 
present Vlasov-MHD model is proper to discuss the solar coronal heating and 
solar wind acceleration by Alfven waves propagating from the photosphere. 
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1 1. Introduction 

2 Large-amplitude Alfvenic fluctuations are ubiquitous in the heliosphere, 

3 especially in the fast solar wind (Bruno and Carbone, 2005). Nonlinear evo- 

4 lution and dissipation of these fluctuations are thought to play important 

5 roles in heating of the solar wind plasmas (Suzuki and Inutsuka, 2006; Wu 
e and Yoon,2007; Araneda et. al, 2008; Valentini et. ai, 2008) and gener- 
7 ation of the localized structures (Tsurutani et. al, 2005; Vasquez et. ai, 
a 2007; Lin et. al, 2009). Since these fluctuations are typically robust for lin- 

9 ear ion-cyclotron damping due to their small wave frequencies and for linear 

10 Landau damping due to their small propagation angle relative to the back- 

11 ground magnetic field, wave-wave interactions (parametric instabilities) are 

12 significant processes of their nonlinear evolution. 
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13 The quasi-parallel propagating Alfvenic fluctuations can resonate with 

14 both the parallel and obliquely propagating magnetohydrodynamic (MHD) 

15 waves (Mjolhus and Hada, 1990; Champouex et. al., 1999; Nariyuki et. al., 

16 2008). In uniform plasmas with a typical parameter of the solar wind, parallel 

17 propagating Alfven waves dominantly resonate with the parallel propagating 
is Alfven waves and ion acoustic waves. It is important that Alfvenic fluctua- 

19 tions can resonate with the "kinetic" wave modes such as ion acoustic waves. 

20 Namely, even if the fluid approximation is valid for Alfvenic fluctuations 

21 themselves, the ion kinetics should be considered for parametric instabili- 

22 ties (Araneda, 1998; Bugnon et. al, 2004; Nariyuki and Hada, 2006; 2007, 

23 Araneda et. al, 2007). 

24 Furthermore, in the plasmas with inhomoginity and/or unstable velocity 

25 distribution functions such as beam components, the resonance and dissi- 

26 pation processes can be different from those in the homogenous plasmas 

27 (Tsikrauri et. al, 2005; Suzuki and Inutsuka, 2005; 2006; Wang et. al, 

28 2006; Suzuki, 2008; Nariyuki et. al, 2009). The kinetic Alfven waves which 

29 propagate at oblique angles relative to the background magnetic field can be 

30 excited by the proton beams in the solar wind (Daughton and Gary, 1998; 

31 Yin et. al, 2007). We note that such beam-excited waves can preferentially 

32 interact with the finite-amplitude parallel propagating Alfven waves. Actu- 

33 ally, the observational studies imply the importance of the kinetic Alfven 

34 waves on the dissipation of the Alfvenic turbulence (Leamon et. al, 1998; 

35 Hamilton et. al, 2008). 

36 As mentioned above, the nonlinear evolution and dissipation processes 

37 of Alfvenic fluctuations in the solar wind are cross scale coupling processes, 

38 in which the MHD-, ion-, and electron-scale phenomena coexist and mutu- 

39 ally interact. Moreover, it is worth noting that the heating and acceleration 

40 process of plasmas by Alfvenic fluctuations are desired to be slot into the he- 

41 liospheric global simulation models (e.g., Nakamizo et. al, 2009). Thus, the 

42 development of the alternative "kinetic" MHD models, which are the "mid- 
43 del scale" model including some non-MHD effects, is necessary for systematic 

44 understanding on such cross scale coupling processes. 

45 As one of such models, we have recently developed a new Vlasov simu- 

46 lation code named Vlasov-MHD code (1-D in the configuration space, 1-D 

47 in "kinetic " velocity space, and 2-D in "MHD " velocity space), in order 

48 to study basic properties of nonlinear evolution of Alfven waves (Kumashiro 

49 et. al, 2009). The concept of the Vlasov- MHD model is the generalized 

50 model of the so-called Landau fluid model (Passot and Sulem, 2003; Bugnon 
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51 et. al, 2004), which includes the nonlinear wave-particle interactions. The 

52 linear analysis of the Vlasov-MHD model has been carried out (Araneda, 

53 1998; Nariyuki and Hada, 2006; 2007, Araneda et. al., 2007) and has con- 

54 eluded that as far as the amplitude of the parent waves is not so large, the 

55 growth rates of the linear analysis are consistent with those in the numeri- 

56 cal results of the ion hybrid simulation. Kumashiro et. al. (2009) performed 

57 the numerical simulation using the Vlasov-MHD code with the Hall-effect 

58 (Vlasov-Hall-MHD code) and demonstrated that the linear growth of para- 

59 metric instabilities of Alfven waves are almost consistent with those of the 
eo ion hybrid simulation. 

ei In the present paper, we discuss the nonlinear evolution of Alfven waves 

62 using the Vlasov-Hall-MHD code. In section 2, we briefly introduce the 

63 numerical schemes of Vlasov-Hall-MHD code. We present the simulation 

64 results in section 3. Section 4 summarizes the results and briefly discuss the 

65 future issues. 



ee 2. Vlasov-Hall-MHD code 

Assuming weak ion cyclotron damping, we include the kinetic effects 
only along the longitudinal (x) direction. Let the ion distribution function 
F(x,t,v) be separated into the longitudinal and perpendicular directions as 
follows 

F(x, t, v) = f{t } x, v x )g(t, x, v y , v z ), (1) 

e/ the Vlasov - Hall -MHD equations are obtained as follows (Araneda, 1998; 
es Nariyuki and Hada, 2006; 2007) 

df Of 1 d / |&J 2 \ Of 

m=- v *dx- + pdx-{ np + —)dV x ' (2) 

du± _ _ u du±_ 
dt dx p dx ' 

db i d ( . i db±_ \ . A . 

u ± -u x b ± --—\. (4) 



dt dx \ p dx 

69 

70 where p is the plasma density (quasi-neutrality assumed), u x is the longitudi- 

71 nal bulk velocity, b± = b y + ib z and u± = u y + iu z are the complex transverse 

72 magnetic field and bulk velocity, and e x is the longitudinal electric field, 
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respectively. All the normalizations have been made using the background 
constant magnetic field, density, Alfven velocity, and the ion gyro-frequency. 

The total pressure is given as p = p e + pi, where p e and pi are electron 
and ion (proton) pressure, respectively. In the present study, isothermal 
electrons are assumed (p e = T e p) (, namely, the total energy in the system 
is not conserved). It is also assumed that the ion and electron pressures are 
isotropic. 

In this study, we solve the Vlasov equation ([2]) with the time-advance 
algorithm called "splitting method" (Cheng and Knorr, 1976), in which the 
Vlasov equation (|2J) is split into the following two advection equations: 



The splitting scheme is widely used because of its simplicity of the algorithms 
and ease of programming. The time advance of distribution function f(x, v x ) 
is first carried out by shifting the distribution function in the x direction (jHJ) 
with the time step At/2, shifting the distribution function in the v x direction 
([6]) with the time step At and again shifting the distribution function in the x 
direction (jSJ) with the time step At/2. The spatial profiles of plasma density 
p is computed by integrating the distribution function over v x . In parallel 
with solving the Vlasov equation, the transverse momentum equation ([3D and 
the induction equation (J3J) are solved by the rational Runge-Kutta scheme for 
time integration and the spectral method for evaluating spatial derivatives. 
The number of cells is 1024 ~ 2048 in the x direction and is 600 in the v x 
direction over a velocity range from v max = 2.4 to v m in = —2.4. The grid 
spacing is equal to Ax = 0.25, and the time step is equal to At = 0.01. 
The boundary condition is periodic for the configuration space and the free 
boundary for velocity space. We adopt PIC scheme (Positive Interpolation 
for hyperbolic Conservation laws) suggested by Umeda (Umeda, 2008) for 
time advancement of the Vlasov equation. 

To analyze the parametric instability of Alfven waves, we give monochro- 
matic, circular polarized and parallel propagating parent waves as initial 
conditions. The initial Alfven wave is written as b± = &o ex P {—ikox), u± = 
uq exp (— ikox), where bo is the amplitude of parent Alfven waves, uo = 
—bo/v^o (Walen relation), where the phase velocity = H>Ao(> 0), Uq = 
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105 A;q(1 + Uq). We adopt the notation that the positive and negative u cor- 

106 responds to the right hand polarized (RH-) and left hand polarized (LH-) 

107 waves. The plasma density is p — po — 1, u x — u x q = 0, and f — fo 

108 is given by the isotropic Maxwellian distribution function. Superposed with 

109 the parent wave is a small amplitude white noise with < \p no i se 

12 > l/2 = 1Q -5_ 

no In the present paper, the ion hybrid simulations are also performed in a 

in way very similar to that described in Nariyuki et. al. (2007) to compare the 

n2 results with those obtained by the Vlasov-MHD simulation. 

n3 3. Simulation results 

iu We first show the simulation results (Run 1) of the parametric instabil- 

ii5 ity of circularly polarized Alfven waves with k = —0.4049, b = 0.25, and 

lie Pi = 0.07, (3 e = 0.5, which parameters are same as those in Araneda et. 

in al. (2008). With these parameters, the modulational instabilities are dom- 

n8 inant (Araneda et. al, 2008). Figure DJa) shows the snapshot of the ion 

n9 distribution at t — 1000, when the instability is almost saturated. Same 

120 as the ion hybrid simulation in Araneda et. al. (2008), ions trapped by the 

121 nonlinear density fluctuations are observed. Figure. QJb) shows the scatter 

122 plot of the ion hybrid simulation plotted in the x - v x phase space at t = 620, 

123 which corresponds to the Fig. 4(c) in Araneda et. al. (2008). 

124 Figure. [2] shows the wave power history of the wave mode with the max- 

125 imum growth rates in the Vlasiv-MHD simulation and ion hybrid simula- 
te tion, respectovely. The linear growth of the Vlasiv-MHD simulation agree 

127 well with the those in the linear analysis, which maxmum growth rate is 

128 0.0141. Furthermore, since the initial parent Alfven wave is weakly nonlin- 

129 ear (6q = 0.0625), the linear growths of the ion hybrid simulation agree with 

130 the linear growth of Vlasiv-MHD simulation as expected by the results in 

131 the past studies (Bugnon et. al, 2004; Nariyuki and Hada, 2007). The dif- 

132 ference of the wave power at early time is due to the numerical noise of the 

133 ion hybrid simulation, which uses the super-particles. 

134 Figure. [3] shows the time evolution of the magnetic field (b±) power spec- 

135 trum, plotted in the phase space of the wave number (k) and time, in the 

136 Vlasov-MHD simulation and the ion hybrid simulation, respectively. While 

137 the modulational instability (around k ~ —0.2 and —0.6) are very similar, 

138 the daughter waves related to the secondary decay instability (wave excita- 

139 tion at k > 0) are excited at the different wave number (around k ~ 0.1 in 
uo the Vlasov-MHD simulation and k ~ 0.05 in the ion hybrid simulation). It 
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141 is worth noting that while the assumptions in the Vlasov-MHD simulations 

142 are thought to be the main cause of such a difference, numerical noises in 

143 the ion hybrid simulations possibly influence the resulting plasma conditions 

144 at the nonlinear stage, since these also influence the parametric instabilities 
us themselves. Thus, the extended Vlasov-MHD simulation and the full-Vlasov 
U6 simulation are necessary to clarify the cause of such differences. 

147 As a matter of fact, the Vlasov-MHD system corresponds exactly with the 

us full-Vlasov system (the system of ion hybrid simulations) under the "static 

149 approximation" p oc \b\ 2 (Mj0lhus, 1976; Mj0lhus and Wyller, 1988), which is 

wo consistent with the condition that parallel electric field (E x ) is the potential 

151 force and the plasmas are isothermal. We here emphasize that in spite of a 

152 lot of assumptions, as far as the nonlinearlity in the system is not so large, 

153 the present Vlasov-MHD simulation model validly solves time evolution of 

154 the Alfven waves both in the linear and nonlinear stage. 

155 We finally mention about the simulation results (Run 2) of the parametric 

156 instability of circularly polarized Alfven waves with k = 0.196, bo = 0.4, and 

157 fa = 0.01, (3 e = 0.32. The decay instability is dominant with these parame- 

158 ters (see Nariyuki and Hada, 2007). Figure H] shows the wave power history 

159 of the wave mode with the maximum growth rates in the Vlasiv-MHD sim- 
leo ulation (black lines) and ion hybrid simulation (gray lines). In contrast to 
lei the previous run, since the amplitude of the initial parent wave is relatively 

162 large (6q = 0.16), the linear growth rates of the Vlasov-MHD system disagree 

163 with the ion hybrid simulation (Nariyuki and Hada, 2007). At the nonlinear 

164 stage, the backward acceleration is caused by the steepen wave packets of ex- 

165 cited daughter Alfven waves in the Vlasov-MHD simulations (Fig. [5](a)). We 
lee remark that such an acceleration is hardly observed in the present ion hybrid 

167 simulation (Fig. [5](b)) . It is because that since the number of accelerated 

168 particle is very small, the number of particles (1500 per cell in the present 

169 simulation), which is not smaller than that in the past studies, is insufficient 
no to clarify the acceleration. 

ni Actually, as shown in Fig. HJ the results of Vlasov-MHD simulations are 

172 not valid when the amplitude of the parent Alfven waves are relatively large. 

173 On the other hand, we infer that such a backward acceleration can occur 

174 in the full Vlasov system. Of course, it does not ensure the quantitative 

175 validity of the Vlasov-MHD model. We believe that the Vlasov-MHD model 

176 can make the suggestions to the ion hybrid simulation, which usually need 

177 much more computational load than the Vlasov-MHD model due to use of 
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Figure 1: (a)The ion distribution function of the Vlasov-MHD simulation plotted in the x 
- v x phase space at t = 1000 (Run 1). (b)The scatter plotts of the ion hybrid simulation 
plotted in the x - v x phase space at t = 620. 
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Figure 2: Wave power history of the wave mode with the maximum growth rates in the 
Vlasiv-MHD simulation (black lines) and ion hybrid simulation (gray lines) . The solid, dot- 
ted and dashed lines indicate the wave power history of the wave mode with k = —0.1593, 
k = —0.1718, k = —0.1847, respectively. The growth rates (wave number) measured in the 
Vlasiv-MHD simulation are 0.01395 (k = -0.1593), 0.01385 (k = -0.1718), and 0.01388 
(k = —0.1847). Those measured in the ion hybrid simulation are 0.01287 (k = —0.1593), 
0.01352 (k = -0.1718), 0.014 (k = -0.1847), which are very similar to those in Araneda 
et. aZ.(2008). 



ns the large number of particles. 

179 4. Conclusion and Discussion 

180 In the present paper, we discussed the nonlinear evolution of finite ampli- 

181 tude Alfven waves using the Vlasov-MHD simulation code, which is a gener- 

182 alized Landau fluid model including the nonlinear wave-particle interactions. 

183 It was confirmed that while the Landau damping should be evaluated along 

184 perturbed field lines (Finn and Gerwin, 1996), as far as the nonlinearlity of 

185 parent Alfven waves are weak, numerical results of the present Vlasov-MHD 

186 simulation agree well with those of the ion hybrid simulations both in the 

187 linear and nonlinear stage. It is worth noting that while the computational 

188 load of Vlasov-MHD code is much less than that of ion hybrid simulation. 

189 On the other hand, when the nonlinearlity of parent Alfven waves are not 
wo weak, the present Vlasov-MHD simulation is not proper. Furthermore, the 

191 Vlasov-MHD model is also inadequate for the high ion beta plasmas, since 

192 the finite Larmor radius effects are neglected. 

193 We note that in spite of restriction on nonlinearlity and ion kinetics, the 

194 present Vlasov-MHD model can be applied to several heliospheric problems. 
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Figure 3: Time evolution of the magnetic field (b±) power spectrum, plotted in the phase 
space of the wave number (k) and time, in (a)the Vlasov-MHD simulation and (b)the ion 
hybrid simulation. 
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Figure 4: Wave power history of the wave mode with the maximum growth rates in the 
Vlasiv-MHD simulation (black lines) and ion hybrid simulation (gray lines). The solid 
and dotted lines indicate the wave power history of the wave mode with k = —0.098172, 
k = —0.12271, k = —0.1847, respectively. The growth rates (wave number) measured 
in the Vlasiv-MHD simulation are 0.0344 (k = -0.098172) and 0.0339 (k = -0.12271). 
Those measured in the ion hybrid simulation are 0.0262 (k = -0.098172) and - 0.0243 
(k = -0.12271). 

195 One of them is the solar coronal heating / solar wind acceleration by Alfven 
we waves propagating from the photosphere (Suzuki and Inutsuka, 2006). Since 
197 the amplitude of the magnetic fluctuations is not so large and the beta ratio 
we is small (Suzuki and Inutsuka, 2006; Tanaka et. al, 2007), the assumption 

199 in the Vlasov-MHD system is well satisfied. 

200 The characteristics of the observed velocity distribution function near the 

201 sun (e.g.Marsch, 2006) is one of the important constraint on the solar coronal 

202 heating and solar wind acceleration model. Namely, the generation processes 

203 of the heliospheric nonequilibrium plasmas such as ion beam components, 

204 temperature anisotropy, and the relative speed among the ion spices should 

205 be comprehensively examined in the model. We note that such a point of 

206 view is important for the collaborating works among future missions (Bepi- 

207 Colombo, SCOPE/Cross-Scale, Solar-C, Solar Orbiter, and so on), which 

208 achievements will contribute toward the heliospheric science. 
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Figure 5: (a)The ion distribution function of the Vlasov-MHD simulation plotted in the x 
- v x phase space at t = 1000 (Run 2). (b)The scatter plotts of the ion hybrid simulation 
plotted in the x - v x phase space at t — 800. In the present paper, 25, 000 particles (1/200 
of the total particles) are plotted. 
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